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Abstract 

We rigorously prove the existence of directed transport for a certain 
class of ac-driven nonlinear one dimensional systems, namely the genera- 
tion of transport with a preferred direction in the absence of a net driving 
force. 

1 Introduction 

In this paper we study the directed transport in ac-driven nonlinear one di- 
mensional systems, whose dynamics is governed by an evolution equation of the 
following type: 

x + jx + U'{x)^E{t). (1.1) 

This equation describes, for example, the classical motion of a particle in the 
potential U{x) under the influence of friction (7 > 0) and space-homogeneous 
external field E(t). We assume U{x) — U{x + 2Tr) a space-periodic function and 
E{t) — E(t -|- T) a time-periodic function with zero mean. 

The paradigmatic case of (jl.ip is the non-conservative nonlinear pendulum, 
U{x) — — cosx, which has been extensively studied. We recall in particular 
the comprehensive review [7] by Mawhin on the forced (damped or undamped) 
pendulum, where one can find a large number of results on the existence of 
periodic solutions of first and second kind, on the existence of Mather sets, and 
on the chaotic behavior of solutions. We also recall the paper [2] by Casdagli, 
dedicated to the study of attractors of dissipative twist maps, where the stro- 
boscopic map of the (periodically) forced damped nonlinear pendulum is given 
as the main example. 

The transport properties of (jl.ip with 7 > have been considered by several 
authors in the physical literature [3-6] . These numerical and theoretical studies 
show that if suitable space-time symmetries are broken then the system can 
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sustain directed current. That is, an ensemble of trajectories with zero current 
at time t = produces a nonzero finite current as t — > c». 

The first paper on this subject [5] was motivated by previous work on 
stochastic ratchets, where directed transport of Brownian particles in asymmet- 
ric periodic potentials (ratchets) is induced by the application of non thermal 
forces. We refer to [8] for an exhaustive review on Brownian motors. In [5] the 
authors study the effect of finite inertia for the ratcheting mechanism in the 
absence of stochastic forces. By using (jl.ip with E{t) harmonic and U{x) a 
two-harmonic non symmetric potential as a working model, they investigate the 
role of regular and chaotic motion in directed transport; in particular, to what 
extent deterministically induced chaos resembles the role of noise. They showed 
numerically that the system can exhibit current flows, whose direction can be 
controlled by varying the amplitude of the external force. See also [6], where 
the origin of this current reversal is identified as a biforcation from a chaotic to 
a periodic regime. 

The role of space-time symmetries is systematically analyzed in [4]. A nec- 
essary condition for the generation of transport is that some relevant space-time 
symmetries in (jl.ip are broken. In fact, nonzero average currents cannot take 
place in the presence of symmetries which allow to generate out of a given tra- 
jectory another one with reversed velocity. In the dissipative case 7 > this 
happens if the potential U{x) is an even function and the driving force E{t) 
is antisymmetric around T/2. In fact, with these hypothesis, (jl.ip is invariant 
under the symmetry (a;,t) {—x,t + T/2). On the contrary, directed current 
may occur if the above symmetry is violated. This can be realized by means 
of an asymmetric driving force [underdamped asymmetrically tilting ratchet) or 
by considering an asymmetric potential [underdamped rocking ratchet). In fact, 
directed current is numerically observed in [4], where, in the case of small dissi- 
pation, the rectification mechanism is explained in terms of a desymmetrization 
of the basins of attraction of two limit cycles with opposite velocities. 

Purpose of the present paper is to rigorously prove the generation of directed 
current due to broken space-time symmetry in (jl.ip . Due to the complexity of 
the dynamics, it is not conceivable to obtain results in the most general case. 
We instead restrict our analysis to the simplest situation when the changes of 
the external field E{t) in the course of time are extremely slow with respect 
to the characteristic time scales of the internal (non-driven) dynamics. In this 
regime, the system is expected to be well described by the so-called adiabatic 
approximation, in which the current at time t has essentially the same value of 
the steady state current Je corresponding to a static external field E = E{t), 
the time t merely playing the role of a parameter. In this approximation, the 
asymptotic time averaged current due to the periodic driving E(t) is simply 

given by ^ dtJE(t)- 

But the above approximation is valid for single valued response functions 
(i.e. not containing hysteretic loops). In the case of (|l.ip this requirement 
turns out to be satisfied if \E\ > max|?7(2:)|, since for such static fields the 
corresponding autonomous system exhibits a unique globally attractive limit 
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cycle. Then, in the case of slowly varying, time-periodic external field which 
fulfills the above condition for large part of the time (see the next section for 
the precise definition), we obtain upper and lower bound on the time averaged 
current by means of its adiabatic approximation. However, this approximation 
is still too complicated to infer the sign of the current, even in special simple 
cases. We then restrict ourselves to the case of strongly enough driving force, 
which allows us to predict the direction of the current in typical examples of 
both asymmetrically tilting and rocking ratchets. 

The plan of the paper is the following: the notation and results are stated in 
Section ^ the steady state current is computed in Section [31 and the adiabatic 
limit is proved in Section SI Finally, Section [S] is devoted to the asymptotic 
expansion of the current for strong external force. 



2 Notation and statement of the result 

Hereafter in the paper we consider the equation 

x + jx + U'{x) ^ Es{Xt). (2.1) 

The potential U{x) is a regular non-constant periodic function of period 27r and 
we set 

- m := min U'{x), M := max U'{x) (2.2) 

(clearly to, M > 0). The positive parameter A appearing in (|2.ip gives the ratio 
between the time scale of the autonomous non-driven system and the one of the 
driving force Eg (r) , which is chosen to be a mollified version of a time-periodic, 
mean zero, piecewise constant function. More precisely, Es{t) is the periodic 
function of period T such that 

^iEi-E2) + ^iEi+E2)Mr~Ts) if t e [0,Ts + 6], 
Es[t) = { (2.3) 
i(Si - E2) ~\{Ei+ E2)Mr -T + S) if r e (T5 + S,T], 

where < 26 <T — Tg, ^^(t) = (j){T/6) with (j){s) a smooth decreasing function 
such that <j){s) — 1, respectively, (f>{s) — —1, if s < 0, respectively, s > 1. The 
time Tg is chosen so that Eg{-) has zero average. 



dTEg{T)=0. 



Note that this implies 



To = lim Tg = — — — T. (2.4) 

" s^a " E1+E2 ^ ' 
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In order to state precisely our result, first we briefly discuss the qualitative 
behavior of the autonomous equation (static E) , 



x + jx + U'{x) = E. 



(2.5) 



We denote by v the velocity of the system, the pair [x, v) thus giving coordinates 
in the cylindrical phase space x R. More precisely, we consider the develop- 
ment of the cylinder on the {x, v) plane. Equation (|2.5p is then equivalent to 
the first order system 

' (2.6) 



E~ U'{x) - jv. 

Note that the equation of the corresponding integral curves on the cylinder is 



dv 

dx ^ 



E - U'{x) 



(2.7) 



Proposition 2.1 Recalling (|2.2p . the following holds. If E > M, respectively, 
E < —TO, there exists a unique globally attractive limit cycle encircling the 
cylinder, which lies in the v > 0, respectively, v < region. We denote the 
cycle by ve{x), that is ve{x) turns out to be the unique periodic solution to 
(j2.7|) . and let Ve be the corresponding mean velocity, that is 



Ve = 27r 



dx ■ 



1 



VEix) 



(2.8) 



Then, there exists C > such that, if x(t; xq^vq) is the solution to (|2.5p with 
initial conditions {xo,Vq) G 5"'" x R, 



x{t;xo,vo) 



t 



Ve 



< C 



(1 + |fo|)log(e+ |i;o|) 
t 



(2.9) 



An appropriate tool for describing transport is to analyze the statistical 
solutions of the dynamics, that is the collective behavior of an ensemble of 
independent particles, whose evolution is governed by the equation of motion 
(j2.1|) . We denote by X\,s{t; xo,vo) the solution to (12.11) with initial conditions 
X\^siO', xq,vq) — xq, X\^s{0',xo,vo) — vq. Given any initial probability measure 
^{dxQ, dvo) on the cylindrical phase space x R, the current at time t is defined 
by averaging the drift of a single trajectory over the initial conditions: 



JxAt) = 



ti{dxo dvg) 



t 



(2.10) 



We say that the system exhibits directed transport if J\,s{t) is definitively away 
from zero as t ^ oo. We will show that this is the case for suitable choices 
of the potential and of the driving force, provided the parameters A and S are 
sufficiently small. More precisely, we first prove that, under further assumptions 
on the driving force, the adiabatic limit A of the current can be explicitly 
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computed. We then show the appearance of nonzero currents in this Hmit 
for suitable choices of the potential and driving force. By continuity, small 
deviation from the adiabatic limit will change numbers but not the fact of 
nonzero currents. 

Theorem 2.2 Recall the definition (I2.10p and assume that 

El > M, E2 > m, (2.11) 

fi{dxo dvo) (1 + \vo\) log(e + \vo\) < 00. (2.12) 



Then there exists the limit 

J = lim lim inf lim inf Ja 5 (i) = lim limsuplimsup Ja ^(i). (2.13) 

Moreover 

J=M^1±I^. (2.14) 

E\ + E2 

3 The autonomous case 

In this section we prove Proposition 12.11 see also [1, Chapter VII, §3], where 
the particular case of the non-conservative pendulum with constant torque is 
treated. We prove the claim for E > M. The case E < —m can be deduced by 
the previous one by noticing that y(t) :— — a;(t; xq, vq) solves the equation 

y + jy + U'{y) = -E 

with initial conditions (— xq, —vq) and potential U{y) ~ U{—y). 

First of all we observe that the divergence of the vector field on the right 
hand side of (|2.6p is equal to —7 < 0. Therefore, by Green's theorem, the 
system (I2.6P does not have any closed paths not encircling the cylinder and can 
have at most one limit cycle encircling the cylinder. We now prove that such 
a limit cycle does exist, it is globally attractive, and lies entirely on the upper 
half w > of the cylinder. 

Given {xq^vq) G S*^ x R, we use in the sequel the abbreviate notation 

x(t) ^ x{t;xo,vo), v{t) = x{t;xQ,vo). (3.1) 
By Duhamcl's formula, 

v{t) ^ e^-^'vo + f dse-^(*-") [£;-[/' (x(s))], 
Jo 

whence 



E-M / E ~ M\ E + m , 

+ Uo e-^* < v{t) < + vo 

7 V 7 / 7 
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It follows that for any a > the strip 



T — 



E-M 



E + m 



1 1 

is positively invariant. In particular, for a suitable number Cq > 0, 

< C'o(l + |wo|) Vt>0. (3.2) 

Moreover, for each a > there exists a positive number Ci = Ci{a) such that 

{{x{t),v{t))^n^ Vi>Cilog(e+|wo|). (3.3) 

Since E > M there are no singular points, therefore the Poincare-Bendixson 
theorem implies the existence of a limit cycle inside the strip TZa- We have 
already shown that this cycle is unique, and therefore asymptotically stable 
with TZa inside its basin of attraction. By (|3.3p we conclude that the cycle is 
globally attractive. We denote by u_e(x) the corresponding integral curve, which 
is therefore the unique periodic solution to (|2.7p . 

It remains to prove the estimate (|2.9p . We fix a < j^^{E — AI) and consider 
the Poincarc map 



G{v) := $(27r;u). 



where ^{x;v) denotes the solution to (|2.7p with initial condition (f>(0;u) = v. 
Clearly v* :— ve{0) is the unique fixed point of G. Since 



dx dv 



dv 



then 



v) 

dv 



exp 



dy 



E-U'{x) d<^{x;v) 
$(a;; vY 

E - U'{y) 



9$(0;w) 
dv 



$(y;w) 



$(a;; v) 



exp 



dy 



<P{y;v) 



where in the last equality we used again that $(a;; v) solves (|2.7p . Therefore 



G'{v*) 



a$(27r;w) 



dv 



exp 



27r 



dx 



ve{x) 



Since ve{ ) > then < G'{v*) < 1, whence there exist 6*2,/? > such that 



\G''{v)-v*\ < Cse 



Vvela VfceN. 



(3.4) 



Recalling the notation (|3.ip . we denote by io, ^1,^2, • ■ •, the increasing sequence 
of times such that x{tk) = (mod 27r) and v{tk) £ Xa (clearly x{tk) ~ x{tk-i) = 
27r). Given i > let n = n{t) = maxjfc >0:tk<t}. By and dSSl), since 
is contained in the upper half f > of the cylinder and xq G [0, 27r), there 
exists C3 > such that 



\x{to)\ < (73(1 + |t;o|)log(e+ |wo|), to < C3log(e+ |wo| 



(3.5) 
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and 

t-t„<C3, X{t) -X{tn) <2tT. 

Noticing that x{tn) — x(to) — 27m we have 



x{t) 



t 



Ve 



< 



\x(t) - x{tn)\ + \x{to)\ + Vsit ~ t„ + to) 



t 



Ve 
t 



2Tm 



where 



27rn 



{tn — to) 



2tt 



E 

n 

E 

k=l 



{tn — to) 



(tk-tk-l) 



dx 



1 



1 



ve{x) <^{x]v{tk-i)) 



(3.6) 



(3.7) 



Since v{tk-x) — G*^ ^{v{to)), by p.4p and the uniform continuity of $(a;;w) in 
TZa^ we get, for some C4 > 0, 



27rn 



(in — ^o) 



< C4 



(3.8) 



Plugging the estimates (|3.5p . p.6p . and p.Sp into (|3.7p we get (|2.9p for a suitable 
positive constant (7. 



4 The adiabatic limit 

In this section we prove Theorem l2.2l By the same reasoning of the previous sec- 
tion, using Duhamel's formula, we conclude that there exists a positive number 
A (independent of A and 5) such that 

\XxAt\xo,vo)\<A{l^\vo\) Vt>0. (4.1) 
Now we define, for any integer fc > 0, 

4'^ = A-ifcT, = \-\kT + Ts), 

tf =X-\kT + Ts + S), tf ^\-\kT + T ~5), 

and denote by iV = N{t) the largest integer k for which t']^ < t. In the sequel 
we shorthand ^^.^(t; a;o, wq) by X{t). We decompose 

N-l 4 

x(t) = x{t) x(tW) + E E \x{tt'') - x{tf) , 

fc=0 i=l 
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with the convention = Since i - t < X'^T and " = 4^^ - 

t'j^'^ ^ A^M, by (gUD we have that 

\X{t)^X{t^j!,^)\ < A(l + |«o|)A-iT, 

|X(4=^)) - X(f ))| + \Xit^^) - X(4^))| < 2A(1 + \v„\)\-'S. 

On the other hand, Es{Xt) = for t e [4^\4^^] ^^^^ i;5(Ai) = -E2 for 
*e [4'\4^^]- Therefore, since 4'^ ~4'^ =A-ir5, 4*^-4'^ =^"'(r-T5 -25), 
and recalhng the hypothesis (|2.1ip . by (|2.9p and (|4.ip it fohows that 

|X(4'^) - X(4'') - X-^TsVe, \<Ci{l + A + A\vo\) log(e + A + A\vo\), 

\X{ti^'>)-X{t'^^^)~X-\T^Ts-2S)V-E, \ < C2 il + A+A\vo\)log{e+A+A\vo\), 



where Ci, respectively, C2 are equal to the number C = C{E) appearing in 
(|2.9p for E = El, respectively, E = —E2- Collecting all the previous estimates 
we conclude that, for some number C > 0, 



Noticing NiXty^ T^^ as t +00 and recaUing ([^ . by assumption fm\i 
the claims (PT5)) and (PH)) follow. 

5 Currents 

In this section we show how directed current actually occurs for suitable choices 
of the driving force and/or the potential. As discussed in the Introduction, 
a necessary condition for the generation of transport is that the space-time 
symmetry (a;,t) — > (— a;, i + A~^T/2) in equation l|2.ip is broken. This symmetry 
is satisfied if the potential U (x) is an even function and the driving force Es (r) 
is antisymmetric around T/2 (this is obtained, see (|2.3p . by choosing Ei = E2, 
Ts = (T— (5)/2, and0(s) antisymmetric around 1/2). It can be violated by means 
of an asymmetric driving force or by considering an asymmetric potential. Since 
the adiabatic approximation (j2.14p is not explicitly known for generic potentials, 
we cannot establish the direction of the current for any value of the parameters. 
We instead expand this expression for large values of the driving force, in order 
to answer the question in at least some special cases. More precisely, we set 
El = E + A, E2 = E and perform an asymptotic expansion of the current for 
diverging values of E. 

Let us first rewrite the expression of J in a more useful form. Recalling 
${x;v) denotes the solution to (|2.7p with initial condition <I>(0;u) = v, we note 



X{t) N 



t xt 



<C- A + <5+- (l + |t;o|)log(e+|«o|). 
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that ^{x,v) := — solves 

dv -E - U'(x) ~ , , , , 

— = -7 + ^, U{x):=U{~x), 

dx V 

with initial condition 'I'(0, v) = —v. In particular, for E > m the unique periodic 
solution to the above equation is ve{x) :— —V-e{~x). Therefore, denoting by 
Ve the mean velocity on the cycle ve{x), the current reads 

J ^ EVe+a - {E + A)Ve 



2E + A 



Ve+aVe 
2tt{2E + A) 



2tt 



E f^'' , E 



dx — / dx 



ve{x) Jo VE+Aix) 



Now we need an asymptotic expansion of the limit cycle for E large. Recall 
ve{x) is a solution to (|2.7p . so that 

~ve{x)^ ^--/VEix) + E-U'ix). (5.1) 
2 dx 

Integrating along the circle, since ve and U' are periodic, we obtain 

I f^"" , , . E 



— dxvE (x) = — 
27r Jo 7 

For any fixed integer iV > 1 and E large enough, we then look for an expansion 
of the following type. 



/N E ^ Vkjx) Rn,e{x) 
"^(^) = - + 2^ + E^+^ ' ^ 

' k=l 

where Vk{x) and RN,Eix) are periodic functions with zero average and Rn,e{x) 
is uniformly bounded as E diverges. By inserting this expansion in (|5.ip and 
equating to zero the first — 1 powers of 1/E, we get that the set of functions 
{vk{x)} must be solution to the following linear differential system in triangular 
form. 



f v[{x) = -^U'{x), 
v'2{x) = -7^«i(a;), 



k-l 



(5.3) 



v'k+i{x) = --y^Vkix) - j^v'g{x)vk-e{x) fc = 3, . . . , TV - 1. 



Once the functions Vk{x) are computed, the remainder Rn,e{x) is determined 
by imposing that the right hand side of (|5.2p is solution to (|5.ip . By a standard 
application of the Gronwall lemma, it is easy to verify that Rn,e(x) turns out 
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to be uniformly bounded as E diverges, we omit the details. Analogously to 
Ve{x), the function ve{x) admits the expansion 



ve{x) = ^ + Y, 

' k=l 



N , 



Vkix) Rn,e{x) 



where the set of functions {vk{x)} satisfy the same system ()5.3p with U'{x) 
instead of U'(x) in the first equation. We note that this implies 



Vk{x) 



Vk{—x) if k is odd, 
— Wfe(— x) if k is even. 



(5.4) 



In the sequel, to simplify the notation, we shall assume, without loss of 
generality, that U (x) has zero average. By using the above expansions for 
= 1 and recalling vi (x) has zero mean, we have 



27T 



dx 



1 



ve{x) 



1 l"^^ 
— dx 



1 



EJo 1 + ^{v,{x) + ^R,,e{x)) 



27r7 7^ 



dxvi{xY + O 



Analogously, 



27r 



dx 



ve{x) 



E 



E^ 



Now, by (|5.3p and (|5.4p . vi{x) — ~jU{x) and vi(x) = vi(—x), whence 



dx vi{xy 



dx vi (x)^ — 7^ 



dxU{xf 



which is positive provided U{x) is not identically zero. Then 



J = 



2tt{2E + A) 



1 



1 



E-i {E + Ay 



2tt 



dxu{xy + o[ 



which implies that for E sufficiently large the current has the sign of A. 
If instead A = the current reads 



J = 



VeVe 
47r 



dx 



1 



1 



VEix) ve{x) 



By the expansion of ve{x) and ve{x) for A^ = 4 and using that Vk{x),Vk{x) 
have zero mean, again by (|5.4p . we have 



dx 



1 



1 



ve{x) ve{x) _ 



27r 



dx 



1 



1 



ve[-x) 



VEix) 



dx Vi (x) V2 (x) + 



67^ 

E» 



dx vi{xy V2{x) 



473 /■^'^ r 1 M 

dx[v2{x)v3{x)+vi{x)v4{x)\ +0 i-^^ 
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Now, by (jO)) . 



dx V\{x) Viix) = \ dx V2{x) V2{x) = 0, 

1 Jo 

dx Vi{x f V2{x) ^ ^ / dxv2{x)'^ V2ix), 



and, integrating by parts. 



Jo 



da; wi (a;) W4 (x) = ^ / da; (x) W4 (x) = ^ / dxv2{x) v'^{x) 

7- Jo 7 Jo 



7^ Jo 



dxv2{x) [7^W3(x) + 7 (uii;2)'(a;)] 



27r 2^ <.2:iT 

dx V2{x) v^{x) -\ — / dxvi{x)v2{x)v'2{x) 
7 Jo 

27r <.27r 

dx W2 (a;) W3 (a;) 5- / c?xw2(x) -y2(a;)^- 

7^ Jo 



Again by (|5.3p . W2(a;) = 7"^G'(x) where G(x) is the antiderivative of U{x) with 
zero average, whence 



dxG{x)G'{xf + o(^-^ 



The integral on the right hand side is nonzero for a generic not even potentiaL 
Let us determine the direction of the current in two typical examples of ratchet 
potential [8]. 

i) Assume U(x) is a piecewise linear "sawtooth" type potential, that is 



C/(x) = 



bx b 

it X e [0,a\, 

a 2 

b{2n-x) b . ^ . 

L _ - if X e a, 2tt] 

27: — a 2 



where a G (0, 27r), & > 0, and the constant term —b/2 makes U of zero average 
(clearly, the equation of motion has to be thought with some mollified version 
of U). The antiderivative of U{x) with zero average is 



G(x) = 



6x^ bx bnin — a) , 
6(27r - x)2 6(27r - x) 67r(7r - a) 



2(27r - a) 
Then, integrating by parts, 



if X G (a, 27r] 



r2TT r2Tr 



dxG{x)G'{xY = -- / dxG(x)'G"(x) = -/i+/2, 
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where 



tt{tt — a) 
3 J ' 

(27r — x)^ 27r — a; 77(77 — a) 1 
2(27r - a) ~^ 2 3 _ ' 

With the substitution x = a + y(27r — a)/a the integral I2 becomes 
so that 

/ dxG(x) G"(x 2 = ^— L / dx { - ^ \-^] ■ 

Jo ^ ' ^ a'' Jo V2a 2 J a 3 J 

Since the integrand is negative for any x G (0,a), we conclude that the current 
is positive for a > 7r and negative for a < tt. In other words, the current has 
the sign of the flatter slope of U{x). 

ii) Assume U{x) is a two-harmonic potential, that is 

f/(a;) = sina; + /xsin(2a;), < /U < 1, 

whence 

G{x) = — cosx — — cos(2a;). 
By an explicit computation, 

j dxG{x)G'{xf = -—n. 

The current is negative, also in this case it goes into the direction of the steeper 

slope of U{x). 

We remark that in the previous examples the current has the same direction 
of that appearing in the analogous cases of overdamped, adiabatically rocked 
ratchets in the zero temperature limit [8]. 
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2a Jo 



2(27r - a) 



X 

2a ' 

27r 



X 

2 



dx 
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